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The helicity dependent parton distributions describe the number density of partons 
with given longitudinal momentum x and given polarization in a hadron polarized longi- 
tudinally with respect to its motion. After the discovery, more than 70 years ago, that the 
proton is not elementary, the observation of Bjorken scaling in the late 1960s lead to the 
idea of hadrons containing almost pointlike constituents, the partons. Since then, Deep 
Inelastic Scattering (DIS) has played a crucial role in our understanding of hadron struc- 
ture. Through DIS experiments it has been possible to link the partons to the quarks, 
and to unveil the presence of other pointlike constituents, the gluons, which lead into 
a dynamical theory of quarks and gluons - quantum chromodynamics (QCD). Polarized 
DIS, i.e. the collision of a longitudinally polarized lepton beam on a polarized target 
(either longitudinally or transversely polarized), provides a complementary information 
regarding the structure of the nucleon. Whereas ordinary DIS probes simply the number 
density of partons with a fraction x of the momentum of the parent hadron, polarized 
DIS can partly answer the question as to the number density of partons with given x and 
given spin polarization in a hadron of definite polarization, either parallel or transverse 
with respect to the motion of the hadron. In here, the phenomenology associated with 
DIS off longitudinally polarized targets, and the present theoretical understanding of the 
dynamics involved, are described. 
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1 Introduction 

Understanding the spin structure of the proton is one of the most challenging problems in 
nowadays subatomic physics. The story of the proton spin dates back to the measurement, 
in 1933, of its anomalous magnetic moment, Hp ^ 1.79 Bohr magnetons, reveahng that the 
proton has an internal structure. We now understand the proton as a bound state of three 
confined valence quarks of spin 1/2, interacting in a comphcated manner through spin-1 
colored gluons. How is the proton built up from its constituents, quarks and gluons, and 
how is its spin distributed among them, are the questions that polarized DIS experiments 
and their theoretical interpretation aim to answer. DIS has always played a crucial role 
in our understanding of the structure of hadrons. The observation of Bjorken scahng 
(Bjorken and Paschos, 1969 [1]) in the DIS experiments in the late 1960s (Breidenbach 
et ai, 1969 [2]) lead to the idea that hadrons contain almost point like constituents, the 
partons. Later, through DIS experiments, it was possible to hnk the partons to the quarks, 
and to discover the existence of electrically neutral constituents, the gluons, which lead 
into a dynamical theory of quarks and gluons - quantum chromodynamics (QCD)(Fritzsch 
et al, 1973 see also H.D. Politzer, 1974 [4]). As a matter of facts, the study of the 
evolution with the momentum transfer of DIS observables represents probably the most 
direct test of the perturbative aspects of QCD. 

Polarized DIS, involving the colhsion of a longitudinally polarized lepton beam on a po- 
larized target (either longitudinally or transversely polarized) provides a complementary 
and important insight into the structure of the nucleon. Whereas ordinary DIS probes 
simply the number density of partons with a fraction x of the momentum of the parent 
hadron, polarized DIS can partly answer the more sophisticated question as to the number 
density of partons with given x and given hehcity, in a hadron of definite polarization, 
either parallel or transverse to its motion. The polarized DIS cross section is described by 
two spin dependent structure functions, gi and g^- The study of the longitudinal polar- 
ization structure function gi for a long time remained comfortably at the level of partons. 
In 1988, the proton data of the European Muon CoUaboration (EMC)(Ashman et ai, 
1988 [5j), which differed significantly from naive theoretical predictions, were pubhshed. 
Those results were argued to imply that the sum of the spins carried by the quarks in a 
proton was consistent with zero, rather than with 1/2, the non-relativistic quark model 
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Figurc 1: Kincmatics of thc DIS process. The quantitics k = {E,k) and k' = {E',k') 
are the four-momenta of the incoming and outgoing leptons, P is the four-momentum of 
a nucleon with mass M and W is the mass of the recoihng system X. We discuss here 
only the exchange of a virtual photon. The four-momentum transferred to the nucleon is 
q = k — k'. 

value, suggesting a spin crisis in the parton modcl. In thc following, thc phcnomcnology 
associated with DIS off longitudinally polarized targcts, and the present theorctical un- 
derstanding of the dynamics involved, are described. At the end, the interested reader will 
find an extended bibliography of up to date professional reviews, and a series of links to 
web pages, where he can satisfy his quest for knowledge and detail, and find the original 
references. 



2 Polarized Parton Distributions 



High-energy lepton scattering off the nucleon, i.e. the process iN i'X, illustrated in 
Fig.l, is called Deep (-g^ = » M^) Inelastic {W^ » M^) Scattering (DIS). The 
filled circle in this figure represents the internal structure of the nucleon which can be 
expressed in terms of structure functions. 

The cross-section for inclusive unpolarized scattering can be written 



dxdy xy 



xy^F,{x,Q^) + 



M'^x^y'^ 



F2{x,Q' 



(1) 



where Fi and F^ are the so called unpolarized structure functions. The following notation 
has been used, 

Q' P-q V .^. 

y = -FT-r. = -f; (2) 
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where i/ = E 



2q-P 2Miy ' " P-k E 
E' is the energy of the virtual photon in the Laboratory frame and x is 
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known as the Bjorken variable. The DIS regime occurs in the so called Bjorken hmit: 



- = oo 



V 



E-E' ^ oo 



X fixed 



(3) 



The difference of the inclusive cross-sections for the lepton and the target nucleon 
polarized longitudinally, i.e. along or opposite to the direction of the lepton beam, is 
given by 



— )■ 



dx dy dx dy 




y^M^x 



, 2M'x'y , 
gi{x,Q ) — — g2{x,Q ) 



(4) 



where the reversal of the nucleon's spin direction is indicated by the double arrow. The 
functions gi^x^Q"^) and g2{x,Q^) are called spin dependent structure functions. Another 
difference of cross sections, similar to Eq. (¡4¡), can be defined for transversely polarized 
nucleons, yielding a different combination of gi{x, Q^) and g^^x, Q"^)- In experiments with 
both longitudinal and transverse target polarization, both gi and g^ can be therefore 
measured. In Eq. ^ it is anyway easily seen that the contribution of g^^x^Q"^) is 
vanishing in the Bjorken limit, Eq. (|3]), and it is therefore difficult to measure it in actual 
DIS experiments. For this reason, only in recent years it has been possible to gather 
precise information on g^^x^Q"^). In this presentation, only the longitudinal polarization 
structure function gi^x^Q"^) wiU be discussed. 



3 The parton model and QCD 

The parton model was proposed by Feynman (Feynman, 1972 [6J), a few years before 
QCD, to explain Bjorken scaling, i.e., the fact that the structure functions in the Bjorken 
limit depend almost only on x and not on Q^. The partons are pointlike constituents 
of the nucleon, which interact electromagnetically like leptons. The nucleon, in a frame 
where it is moving very fast, could be viewed as a "beam" of collinear partons. The 
partons are characterized as having momentum p = x'P, where P is the momentum of 
the nucleon, and covariant spin vector s. The interaction with the hard photon is then 
visualized as in Fig. 2, in which the lepton-parton scattering is treated analogously to 
elastic lepton-lepton scattering. 5* is the covariant spin vector of the nucleon. Requiring 
the final parton to be on mass shell, i.e. (p -|- g)^ = 0, selects the value x' = x. Thus 
X can be interpreted as the fraction of longitudinal momentum of the target carried by 
the struck parton. The partons were shown to have the quantum numbers of the current 
quarks, which were found to carry only a part of the nucleon momentum. The missing 
part was ascribed to a neutral vector particle, the gluon, which only carries color degrees 
of freedom, and therefore does not couple to leptons. 

For unpolarized DIS, one finds the scaling result expressed in terms of the number 
density of quarks, q{x), and antiquarks, q{x) 

^i(^) = ^Ee?[^^(^) + ^i(^)]' (5) 
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Figure 2: Parton model description of DIS. 

where the sum is over flavors j, ej is the charge of the quarks, and the Callan-Gross 
relation 

F2(x) =2xFi(x), (6) 

imphes that the charged partons are spin 1/2 particles. 
For longitudinally polarized DIS one obtains 

9ii^) = ^ E e| [Ag,(x) + Aq,ix)] , (7) 

where 

Aq{x) = q^{x) - q^{x) , (8) 

are the longitudinal polarization functions, or helicity dependent parton distributions, 
with g_j.(.¡_) (x) representing the number densities of quarks whose spin orientation is par- 
allel (antiparallel) to the longitudinal spin direction of the proton (see Fig. 3). In terms 
of these, the unpolarized parton density is 

q{x) = q^{x) + q^{x) . (9) 



P 




Figure 3: Visualization of the longitudinally polarized parton density Aq{x). The upper 
arrows show the spin direction. 

The parton model is an intuitive description, much like the impulse approximation in 
nuclear physics, and it appeared long before QCD. Once QCD is accepted as the theory of 
strong interactions, with quark and gluon as the fundamental fields, there are interaction 
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Figure 4: QCD generalization of the parton model. In the factorized scheme, H represents 
the hard part, while $ the soft part. 

dependent modifications of the simple parton model resuhs for DIS. The main impact 
of the QCD interactions is: i) to introduce calculable logarithmic dependence in the 
parton densities; ii) to generate a contribution of the gluons to the structure functions, 
and, in particular to gi, arising from the polarization of the gluons in the nucleon. 

Unfortunately, these correction terms are infinite. The infinity is caused by collinear 
divergences which occur because of the masslessness of the quarks and which are removed 
through a mechanism called factorization of collinear divergences. This mechanism allows 
the reaction to be written as a product of a hard and a soft part (see Fig. 4), and the 
infinity is absorbed in the latter, which cannot be calculated in QCD. For this part either 
models of hadrons or parametrizations of data are used. After an elaborate calculation, 
gi is found to depend also on Q^: 

9iix,Q') = l H e¡{Aq{x,Q') + Aq{x,Q') 

flavors 

+ /' -{^C,{x/y) [Aq{y, Q') + Aq{y, Q')] 

Zn Jx y 

+ ACG{x/y)AG{y,Q')}} (10) 

where «^(Q^) is the QCD running coupling constant, while ACg and ACq are the so 
called Wilson coefficients, which can be evaluated in pQCD. One can notice that the 
first line in the above equation represents the result obtained in QCD at leading order, 
while the remaining part is given by the higher-order contributions. Both the helicity 
dependent parton distributions and the Wilson coefficients depend on the factorization 
and renormalization schemes used. 
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4 The proton spin problem 



Let us return to the parton model and write the expression of gi, for protons (p) and 
neutrons (n), respectively, in terms of hnear combinations of the hehcity dependent parton 
distributions associated with the quark flavors u, d, s: 



±^Ag3(x) + ^Ag8(x) + AE(x) 



where 



Agg = ( An + Am) - (Ad + Ad) , 
Agg = (Ai¿ + Aü) + (Ad + Ad) - 2(As + As) , 
AS = (Am + Am) + (Ad + Ad) + (As + As) . 



(11) 

(12) 
(13) 
(14) 



These quark densities transform respectively as the third component of an isotopic spin 
triplet, the eighth component of an SU{?>)p octet, and a flavor singlet. 
Taking the flrst moment of Eq. f lTT|) yields 



pp{n) 
J- 1 



/•1 1 

/ gi{x)dx = - 
Jo 9 



^3 3 



19\ + 9a 



(15) 



where 



9l 



9\ 



9\ 



dx Ag^^x) 
dx Aqs{x) 



dx AIl(x) 



(16) 



The octet of currents associated to the distributions Eqs. (ÍT2l) - (ÍT^ is precisely the one 
that controls the weak decays of the neutron and of the octet hyperons, which imphes 
that the values of g\ and g\ are known from /3-decay measurements: 



^3=^^ = 12670 + 0.0035, = 0.585 ± 0.025 . (17) 

Hence, according to Eq. ( ÍTSl) . a measurement of Fi can be considered as giving the value 
of the flavor singlet gj^. 

In the parton model, two remarkable sum rules can be obtained. The Bjorken sum rule 
(Bjorken, 1970 [7J), initiaUy derived using only the parton model and isospin invariance 
(see Eq. (ÍT5|) ). reads, to leading order in QCD: 

í'dx[gUx,Q')-g-{x,Q')] = ^. (18) 
Jo 6 



6 



QCD corrections to this result have been calculated up to third order in the couphng 
constant (Larin et al, 1997 [8]) leading to 

£ dx [gl{x, Q') - Q')] = ^ |l - ^ - 3.583 (^)' - 20.215 (^)'| , (19) 

a result which has been confirmed experimentally. 
Besides, in the parton model 

í?o =(78 +3(As + As), (20) 

and, if one neglects the contribution from the strange quarks in Eq. fl2Dl) . i.e., from 
As + As, the EUis and Jaffe sum rule (Elhs and Jaffe, 1974 [9]) is obtained : 

~ ^ 0.6 . (21) 

In 1987, the EMC Cohaboration performed a measurement of and Ti and, using the 
known values of g\ and g% in Eq. (fTSl) . obtained 

9a ^ 0, (22) 

a value in contradiction with the Elhs-Jaffe sum rule. This experimental result has other 
dramatic imphcations. Consider the physical significance of AZl(a;). Since g_>.(^)(x) count 
the number of quarks of momentum fraction x with spin component ±i along the direction 
of motion of the proton (say the z-direction), the total contribution to Jz coming from 
the spin of a given fiavor quark is 




It foUows that 

^ 2(5f-^^) , (24) 

where is the contribution to Jz from the spin of all quarks and antiquarks. The 

connection between g^ and (5'|"'"''^*) , and the possible implications of the result (5*9"'^^''«*) 
smaller than 1/2, were discussed in a paper by Sehgal [TÜ] . 

In a non-relativistic constituent model one would naively expect all of the proton spin 
to be carried by the spin of its quarks. In a relativistic model one expects 2(S'|""^'^*) ^ 0.6, 
due to the loss of normalization of the upper components, in agreement with the Ellis- 
Jaffe sum rule but far from the EMC result for g% This contradiction was labeled as 
the proton "spin crisis", and it was the beginning of an impressive experimental and 
theoretical activity, which is still going on. From the theoretical side, a careful reanalysis 
both of non relativistic and relativistic models has lead to results in closer agreement 
to the EMC value, as we will show. Another argument has been the observation that. 
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even in the Bjorken limit, the gluonic version of the anomalous triangle diagram leads to 
gluonic contribution to the first moment of gi (Altarelh and Ross, 1988 [TT]): 

rf"°"'(g2) = _1 AG(g2) . (25) 

This result is of fundamental importance since it imphes that the simple parton model 
formula for g'^ is incomplete. Instead, one has that the quantity measured by the EMC 
collaboration is 

fA = 9'A-^^^G, (26) 

aUowing for the EJ sum rule to be fulfiUed for a large enough value of AG. However, 
this resuh depends on the factorization scheme utihzed. In schemes where g\ has the 
meaning of a spin, the small measured g\ does not necessarily imply that the physically 
meaningful g^ is small. This observation was initially presented as a possible resolution of 
the spin crisis, but it is now clear that it is not sufñcient. As a matter of fact, to explain 
the measured value of g% it should be AG ~ 1.7, while the observed values are much too 
small (|AG| ?a 0.29) to resolve the spin crisis, and different analyses give even different 
signs. A way out of this problem is to consider that the partons possess orbital angular 
momentum. 

From the experimental side, after the publication of the EMC data, an impressive 
program started in several laboratories to extend the kinematical range of the EMC 
experiment, to reduce the systematic errors in the measurement of F^, and to get the 
neutron information, necessary to test the fundamental Bjorken Sum Rule, Eq. f lT8|) . In 
this way, it was possible to realize if the spin crisis could be ascribed to a problem of the 
parton model and, in turn, of the underlying theory, QCD. The neutron measurement is 
really difñcult due to the lack of pure neutron targets, so that nuclear targets have to be 
used and nuclear structure effects have to be carefully taken into account. Despite of these 
difñculties, a reasonable amount of precise data is nowadays available for the neutron, as 
it will be seen in the following. 

5 Models of hadron structure 

In this section, it will be shown that properly built models can help in clarifying the origin 
of the so-called spin crisis. QCD is a theory of quarks (antiquarks) and gluons, as has been 
shown in the asymptotic regime, where the interaction can be treated perturbatively. At 
low energies, the idea that baryons are made up of three constituent quarks and mesons of 
a constituent quark-antiquark pair, the naive quark model scenario, accounts for a large 
number of experimental facts. The quest for a relation between the "current" quarks of 
the theory and the constituent quarks of the model has an old history and this search 
has been the leitmotiv of a considerable research effort. The fundamental problem one 
would like to understand is how confinement, i.e. the apparent absence of color charges 
and dynamics in hadron physics, is realized. 
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Figure 5: The structure function F^^x^Q'^) obtained by NLO-evolution to = 10 
GeV^ in a model with point-hke quarks (dashed) and in the convolution approach with 
composite quarks (fuU) using the wave functions of a semirelativistic model. See Ref [12] 
for details. 



Detailed quark models of hadron structure based on the constituent quark concept have 
been defined in order to explain low energy properties. In order to reach the high energies 
of the DIS regime, the central assumption is the existence of a scale, /ig, where the short 
range (perturbative) part of the interaction is negligible, and therefore the glue and sea 
are suppressed, and the long range (confining) part of the interaction produces a proton 
composed of three (valence) quarks only. One then ascribes the quark model calculations 
of matrix elements to that hadronic scale //q- For larger their Wilson coefScients will 
give the evolution as dictated by pQCD. In this way quark models, summarizing a great 
deal of hadronic properties, may substitute " ad hoc" low-energy parameterizations. The 
procedure describes successfully the gross features of the DIS results. 

In order to produce a more quantitative description of the data, different mechanisms 
have been proposed: sea gluons, sea quarks and antiquarks, meson clouds, relativistic 
effects, etc... Some of these mechanisms appear naturally if one endows the constituent 
quarks with structure. In this scenario, the constituent quarks are themselves complex 
objects whose structure functions are described by a set of functions that specify the 
number of point-Iike partons of type b, which are present in the constituents of type a 
with fraction x of its total momentum. In general a and b specify all the relevant quantum 
numbers of the partons, i.e., color, flavor and spin (see Altarelli et ai, 1974 [TB], where 
this scenario has been flrstly addressed. See also Hwa, 1980 [14j). 

Results are shown in Figs. |5]and[6l Fig. |5]refers to the unpolarized case. The structure 
function F^^x^Q"^), calculated evolving the parton distributions obtained at low energy 
Uo and do within a semirelativistic model, describes successfully the data. Moreover, the 
agreement becomes impressive when compared with the analogous calculation with non- 
composite constituents. A qualitatively similar agreement is obtained also in the polarized 
case, as it is shown in Figs. 6 and 7. It should be noticed that in this framework the spin 
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Figure 6: Left (Right): xgi{x, Q"^) for the proton (neutron) evolved at NLO to = 10 (5) 
GeV^, for two different convolution models (full curves). The proton (neutron) data are 
from the EMC and SMC collaborations (from the E154 colaboration at SLAC). See Ref 
[15] for details. 



crisis, as initially presented, does not arise: the constituent quarks, the effective degrees 
of freedom of the model, can carry most of the proton spin, even if the measured g\, a 
current quark DIS observable, is smalL In this scenario, there is no need of large orbital 
angular momentum contributions to explain the data, once the proper structure of the 
constituent quarks has been taken into account. 

In general, the non relativistic constituent quark models describe the nucleon by S- 
wave functions, or with small D-wave admixtures, once the hyperfine interaction is taken 
into account. Relativistic constituent quark models have a non-zero quark orbital angular 
momentum from the beginning and reduce the valence quark spin contributions to the 
nucleon spin from 1 to about 0.6. However, to reach the scale of the data, perturbative 
QCD evolution has to be used and the contribution of quarks and gluons to the orbital 
angular momentum may be of different size, and relevant at the experimental scale. How- 
ever, the perturbatively generated gluon orbital angular momentum is cancelled by the 
gluon helicity contribution, as it can be shown by general arguments. 

6 Experimental determination of the longitudinal po- 
larization functions 

A vast amount of data on the inclusive spin structure function gi{x,Q^) has been accu- 
mulated by lepton scattering experiments at SLAC, CERN, DESY and Jefferson Lab (see 
Fig. 7). 

Inclusive data from proton, neutron, deuteron and perturbative QCD analysis are used 
to extract the contributions from quark, antiquark and gluon densities with a precision 
that depends very strongly on that of the data. One approach to gather additional 
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Figure 7: World data on the polarized structure function xgi (x) for the proton, deuteron 
and neutron in the DIS region {W > 2 GeV), taken by different experiments, at several 
different values of Q^, as compiled in Ref. [I6¡ . 



information has been to use semi-inclusive lepton scattering (SIDIS), where, in addition 
to the scattered lepton, one detects a leading hadron (typically a pion or kaon) in the 
final state. This approach was pioneered by the SMC collaboration and the most detailed 
data set stems from the HERMES and COMPASS coUaborations (see Fig. 8). However, 
with SIDIS alone it is difñcult to extract the sea and glue polarization. An alternate 
route has been employed in the polarized proton collision program at BNL, using the 
Relativistic Heavy lon Collider (RHIC). One selects final state signatures like hadrons, jets 
or even direct photons with large transverse momentum that indicate an underlying hard 
interaction between two constituents from the two colliding protons. Recent analyses take 
into account all the data from inclusive, semi-inclusive polarized deep inelastic scattering 
and polarized pp scattering at RHIC, and perform all theoretical calculations at next- 
to-leading order of perturbative QCD to maximally constrain the extracted distributions 
(De Florian et al., [17] ), opening the door to obtaining a better and more reliable picture 
of the spin structure of the nucleon. 
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Figure 8: World data on helicity dependent parton distributions Aq{x) extracted from 
semi-inclusive DIS data, as compiled in Ref. [Tü] . 



7 Summary 

After 30 years of dedicated experiments, a rather detailed picture of the nucleon spin 
structure is arising. It appears that quark helicities do contribute a significant fraction 
of the (longitudinal) spin of the nucleon. Most importantly, the fundamental Bjorken 
sum rule and its pQCD evolution are consistent with the data. None of the existing 
experiments show any features that would contradict perturbative QCD in its realm of 
applicabihty. 

Most pQCD-based analyses of the data agree fairly well on the contribution of various 
quark fiavors to the proton spin, although there is still some controversy on the role played 
by strange sea quarks. We also do not know yet whether polarized up and down anti- 
quark densities show the same difference as the unpolarized ones. The emergent picture is 
one where valence quarks {Aqv = Aq — Aq) carry roughly the expected fraction (~ 60%) 
of the nucleon spin, while the (on average) negative hehcity of sea quarks reduces this 
to about 30-35 at Jefferson Lab wiU extend our knowledge of polarized quark densities 
out to X > 0.8 and will decisively test predictions from pQCD and QCD-inspired models. 
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Additional information on individual quark and antiquark flavor contributions to the 
nucleon spin will come from future experiments at RHIC (in particular from direct 
production at 500 GeV center-of-mass energy) and the new FAIR facility in Darmstadt, 
Germany as well as semi-inclusive measurements with COMPASS and at Jefferson Lab. 

In addition to logarithmic violations of scaling, expected from pQCD, the data show 
some evidence for the so called higher twist contributions, i.e., non-scaling contributions 
to the structure function gi, at intermediate values. Future measurements at Jefferson 
Lab wiU improve our knowledge of higher-twist matrix elements, and of the second spin 
structure function, g^, and its moments. 

At the lower end of the scale, at the real photon point, the Gerasimov-Drell-Hearn 
(GDH) sum rule (Gerasimov, 1966 [18]; Drell and Hearn, 1966 [19]) has been well con- 
firmed. The transition from this point to DIS, through the region of intermediate Q^, is 
under investigation, in particular at JLab. 

One of the most important remaining open questions is where is the nucleon spin that 
is not carried by quark helicities. We do not know the precise magnitude and shape of 
the gluon contribution, and it could stiU be an important fraction of the total. Improved 
statistics from the direct measurements at RHIC and COMPASS wiU help dramatically. 

The contribution from quark angular momentum is also an important ingredient in 
the total spin balance of the nucleon. While a direct measurement is not available, 
one can learn much about the transverse distribution and motion of quarks from semi- 
inclusive measurements of single spin asymmetries. This is a fairly new fleld, related to 
the transversity observable, with a very rich potential and a rapidly growing body of 
experimental data, but lies outside the scope of this article. 

Another way to describe the nucleon spin is the study of Generalized Parton Distri- 
butions (GPDs), in particular Deeply Virtual Compton Scattering (DVCS). Moments of 
certain combinations of GPDs can be related to the total angular momentum (spin and 
orbital) carried by various quark flavors, as expressed in Ji's sum rule (Ji, 1997 [2Ü]). 
Data taken at Jefferson Lab are being analyzed and future experiments at COMPASS 
wiU clarify the situation enormously. 

In summary, the description of the nucleon spin keeps being a fundamental problem of 
hadron structure (see Refs. [211 [221 [231 [21 ESl ESl EZ] and the links in Ref [2B] for further 
reading). An extensive and rich experimental program has been undertaken at existing 
facilities like CERN (COMPASS), RHIC and Jefferson Lab and much theoretical effort is 
being devoted to a fuU understanding of the nucleon spin structure. 
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